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ABSTRACT 

Massive star formation occurs in Giant Molecular Clouds (GMCs); an understanding of the evolution 
of GMCs is a prerequisite to develop theories of star formation and galaxy evolution. We report the 
highest-fidelity observations of the grand-design spiral galaxy M51 in carbon monoxide (CO) emission, 
revealing the evolution of GMCs vis-a-vis the large-scale galactic structure and dynamics. The most 
massive GMCs (Giant Molecular Associations - GMAs) are first assembled and then broken up as the 
gas flow through the spiral arms. The GMAs and their H2 molecules are not fully dissociated into 
atomic gas as predicted in stellar feedback scenarios, but are fragmented into smaller GMCs upon 
leaving the spiral arms. The remnants of GMAs are detected as the chains of GMCs that emerge 
from the spiral arms into interarm regions. The kinematic shear within the spiral arms is sufficient to 
unbind the GMAs against self-gravity. We conclude that the evolution of GMCs is driven by large- 
scale galactic dynamics -their coagulation into GMAs is due to spiral arm streaming motions upon 
entering the arms, followed by fragmentation due to shear as they leave the arms on the downstream 
side. In M51, the majority of the gas remains molecular from arm entry through the inter-arm region 
and into the next spiral arm passage. 

Subject headings: galaxies individual (NGC5194, M51) — ISM: clouds — ISM: evolution 



1. INTRODUCTION 

Despite numerous studies of molecular gas i n the 
Milky Way an d galaxies (|Scoville fe Sanders! 119871 : 
IBlitz et ail 120071 ). the processes affecting evolution of 
GMCs have remained poorly understood. In fact, 
uncertainty remains as to whether GMCs are stable 
structures th at survive over a galac t ic rotation period 
(> 1 s yrs; IScoville fc Hershl 119791 IScoville fc WilsorJ 
2004) or are transient structures, destroyed immedi- 
ately after formation by violent feedback f rom young 
stars [a few x 10 7 yrs; IBlitz fc Shul 119801: lElmegreenl 
20071 or even shorter ~ 10 6 yrs: lElmegreer] 120001: 



Hartmann. Ballesteros-Paredes fc Berginll2001j . 

The internal structure of GMCs is generally believed 
to be determined by an approximate balance of self- 
gravity and internal turbulent or magneto-turbulent 
pressures (exceeding the e xternal ISM pressures by 
two orders of magnitude; iMverd 119781) . Therefore, 
GMCs are likely self-g ravitationally bound, dis tinct, 
and long-lived objects (jScoville fc Sandersl Il987t ). A 
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difficulty, however, arises in maintaining the internal 
turbulence over the long lifetimes, since the turbu- 
lence should dissipate rapidly within a cloud crossing 
time (~a few xl0 6 yrs). The energy must be con- 
tinuously resupplied if their lifeti mes are longer (e.g. 
iHeitsch. Mac Low fc Klessenl 120011 ) . The energy source 
is still unknown although supernovae and galactic 
rotation have been suggested (|Mac Low fc Klessenll2004l : 
iWada. Meurer fc Normaij2002tlKoda et al.ll2006fl . Mod- 
els for rapid GMC formation and disruption have been 
proposed to avoid the need to re-energize the turbul ence 
( Balleste ros-Paredes. Hartmann. Vazquez-Semademl 
I1999D . The absence of older > 5 Myr stellar populations 
within GMCs has been cited as evidence of short GMC 
lifetim es (jHartmann. Ballesteros-Paredes fc Berginl 
120011) . However, more re cent obs e rvations do find some 
older star populations ( iJeffriesI 120071 : iGandolfi et al.1 
l2008t lOhveira et al.ll2009h . presumably indicating longer 
lifetimes. Moreover, at some point the older stars 
must drift free of the GMCs since they are no longer 
subject to the hydrodynamic forces of the ISM gas. The 
collapse timescale of star-forming cores via ambipolar 
diffusion is substantial, and may support the longevity 
(jTassis fc Mouschoviasll200l . 

This controversy carries over when considering the 
galactic spatial distribution of GMCs, since early molec- 
ular line surveys of the Milky Way linke d the GMC 
evolut ion to the Galactic disk dynamics. iCohen et al.l 
(fl980h argued that GMCs were confined largely to the 
spiral arms with very few seen in the inter-arm re- 
gions, implying that GMCs must be short-lived with 
a lifetime similar to the arm-crossing timescale (a few 
x 10 7 vr: iDame. Hartmann fc ThaddeusI 120011) . How- 
ever, [Sandere^coyin^^^olomoiil W8l|) claimed to find 
many GMCs in the inter-arm regions, suggesting that the 
GMCs must be quasi-permanent structures, surviving 
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> 10 8 yrs, i.e. a substantial fraction of a galactic rotation 
period. A recent 13 CO survey sugge sts that GMCs in 
the inter-arm region s are less-massive (|Koda et al.l l2006; 
Uackson et al.l [2006h . accounting for the earlier discrep- 
ancies if the Cohen et al. survey missed the less-massive 
inter-arm GMCs. In any case, all Galactic surveys have 
the intrinsic limitation that the velocity dispersion and 
streaming motions of the clouds can blur out the spiral 
arm/inter-arm definition when the disk is viewed edge- 
on. 

High-resolution molecular gas observations of external 
face-on spiral galaxies are therefore essential for the full 
census of GMC population over galactic disks. How- 
ever, prior interferometers, which are required for such 
high-resolution imaging, had only a small number of tele- 
scopes, and thus were severely limited by low image- 
fidelity. Indeed, the high side-lobes of bright spiral arms 
due to poor itu-coverage have often led to false struc- 
tures in the inte r -arm regions (IRand fc Kulkarnil[l990l ; 
lAalto et al.lfl99l IHelfer et al.H2003n . ~ 

2. OBSERVATIONS AND DATA REDUCTION 

High-fidelity imaging of nearby galaxies at millimetre- 
wavelengths has now become feasible with the Combined 
Array for Research in Millimeter Astronomy (CARMA). 
CARMA is a new interferometer, combining the six IO- 
meter antennas of the Owens Valley Radio Observatory 
(OVRO) millimeter interferometer and the nine 6-meter 
antennas of the Berkeley-Illinois-Maryland Association 
(BIMA) interferometer. The increase to 105 baselines 
(from 15 and 45 respectively) enables the highest fi- 
delity imaging ever achieved at millimeter wavelengths. 
The entire optical disk of the Whirlpool galaxy M51 
(6.0' x 8.4') was mosaiced in 151 pointings with Nyquist 
sampling of the 10m antenna beam (FWHM of 1 arcmin 
for the 115GHz CO J=l-0 line). The data were reduced 
and calibrated using the Multichannel Image Reconstruc- 
tion, Imag e Analysis, and Display (M IRIAD) software 
package (|Sault. Teuben fc Wright]ll995h . 

We also obtained total power and short spacing data 
with the 25-Beam Array Receiver System (BEARS) 
on the Nobeyama Radio Observatory 45m telescope 
(NR045, FWHM = 15") Usin g the On-The-Fly observ- 
ing mode (Sawada et al.l l2008h . the data were oversam- 
pled on a 5" lattice and then re-gridded with a spheroidal 
smoothing function, resulting in a final resolution of 22". 
We used the NOSTAR data reduction package developed 
at the Nobeyama observatory. We constructed visibili- 
ties by de-convolving the NR045 maps with the beam 
function (i.e. a convolution of the 15" Gaussian and 
spheroidal function), and Fourier-transforming them to 
the uu-space. We combined the CARMA and NR045 
data in Fourier space, inverted the uv data using theo- 
retical noise and uniform weighting, and CLEANed the 
maps to yield a three-dimensional image cube (Right As- 
cension, Declination, and LSR Doppler velocity). 

The combined data have an RMS sensitivity of 40 
mjy/beam in 5.1 km s _1 wide channels, corresponding 
to 1 x 1O 5 M at the distance of 8.2 Mpc (adopting a 
CO-to-H 2 conversion factor of Xqo — 2 x 10 20 cm~ 2 
[K km s" 1 ]- 1 ). Typical GMCs in the Milky Way 
(i.e. 4 x 10 5 M(7) in m ass and 40 pc in diameter, 
IScoville fc Sanders! fl987f) are therefore detected at 4cr 
significance. Our angular resolution of 4" corresponds to 



160 pc, which is high enough to isolate (but not resolve) 
the GMCs, given that the typi cal separation of Galactic 
GMCs is a few 10 pc to kpc (IScoville fc Sanders! [19871 : 
IKoda et al.l 120061 ). The combination of spatial resolu- 
tion, sensitivity, and image -fidelity differentiates our 
study fr om previous work dVogel. Kulkarni fc Scovillel 
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Heifer et al. I 12003ft . and enables a reliable census of 



GMCs in M51. 

Figure QJl shows the CO map integrated over all ve- 
locities. The new image shows the full distribution 
of molecular gas over the entire optical disk of M51 
(14.3 x 20.0 kpc 2 ), including both the prominent spi- 
ral arms and inter-arm regio ns. The bright inner arms 
have been previously im aged (jVogel. Kulkarni fc Scovillel 
fl988t lAalto eTall fl999h . but the new data extend these 
arms over the full disk and most importantly, yield sig- 
nificant detection of inter-arm GMCs for the first time. 
Figure [2^ shows the distribut ion of discrete GMCs mea- 
sured using CLUMPFIND (|Williams. de Geus fc Blit3 
Il994f >. clearly indicating many GMCs with mass exceed- 
ing 4 x 10 5 Mq in the inter-arm regions. Lower mass 
GMCs would not be detected by the cloud finding algo- 
rithm and only 36% of the inter-arm CO emission is seen 
in the detected discrete clouds shown in Figure 2. 

3. DISCUSSION 

Two possibilities for GMC formation and lifetime can 
clearly be distinguished from the observed GMC distri- 
bution relative to the spiral arms. One possibility is that 
they are transient, short-lived structures - formed lo- 
cally at convergence locations of the galactic hydrody- 
namic flows and destroyed quickly by disruptive feed- 
back from star formation within the GMCs. Alterna- 
tively, if the GMCs are long-lived, lasting through the 
inter-arm crossing period, their presence in the inter-arm 
regions is naturally explained. A simple calculation can 
be done to rule out the formation of abundant GMCs 
from the ambient gas with a low average-density in the 
inter-arm regions. Assuming a very ideal spherical gas 
accretion of the converging velocity v into the volume 
with the radius r, the mass accumulates within the time 
t is M = 4irr 2 vmnnt, where n is the number density of 
ambient gas and mn is the mass of hydrogen atom. Us- 
ing a GMC radius r ~ 20 pc, typical velocity dispersion 
in galactic disk v ~ 10 km s _1 , and average gas density 
n ~ 1 — 5 cm" 3 (inferred for M51 in the areas without 
GMCs), it takes ~ 10 8 yrs to accumulate 4 x 1O 5 M . 
This mass-flow argument is valid even when the ambient 
gas is not exactly diffuse but consists of smaller clouds if 
their distribution is roughly uniform. Therefore, the ma- 
jority of the GMCs in the inter-arm regions cannot have 
formed there locally on the required short timescales (i.e. 
inter-arm crossing timescale ~ 10 8 yrs). 

Figure [5^ also reveals that the molecular gas prop- 
erties, specifically their masses, are significantly depen- 
dent on galactic environment. The most massive struc- 
tures with 10 7_8 Mq, referred to as GMAs, are found 
only in spiral arms, not in inter-arm regions. The im- 
proved image fidelity was necessary to avoid misidenti- 
fication of even such massive GMAs, especially in the 
inter-arm regions. GMAs must therefore form and dis- 
rupt while crossing the spiral arms (with a timescale 



Dynamically Driven ISM Evolution in M51 



3 



typically ~ 2-5 x 10 7 yrs). Their formation within the 
arms is aided by the spiral streaming which causes de- 
flection and convergence of the galactic flow streamlines 
in the arms. Although stellar feedback (e.g. photo- 
dissociation by OB star ultraviolet radiation and super- 
nova explosio n s) is often invoked for G MC destruction 
(|Larsonlll987t I Williams fc McKed fl997l) . these mecha- 
nisms are not a likely cause of significant GMA dis- 
sipation given their mas ses (10 7_8 M Q , see Figure [2^; 
iWilriams fc McKedfl997t ). More telling is the high mass 
fraction of H2 in the inter-arm regions (Figure Eb,c). 
Comparing the mean molecular g as surface density with 
that of the HI (|Braun et al.ll2007h . we estimate that 70- 
80% of the gas remains H2 within the major part of the 
disk (~ 12 kpc; Figure [2):). Thus, GMAs are not signif- 
icantly dissociated into atomic or ionized gas. We con- 
clude that the GMAs must be fragmented into the less- 
massive GMCs and the most massive of these are seen as 
discrete clouds in the inter-arm regions in Figure^. Fig- 
ure^ includes only the GMCs down to the 4cr level; the 
sum of their emission is 36% of the total emission of the 
entire disk. The remainder of molecular emission is pre- 
sumably in less-massive GMCs, not identified as discrete 
clouds at the current resolution and sensitivity. Their 
distribution must be roughly uniform in the inter-arm 
regions since they are resolved out at 4" resolution. 

It is unlikely that a significant portion of the remain- 
ing emission arises in diffuse molecular gas at low den- 
sities. First, the critical density for collisional-excitation 
of CO(J=1-0) transition is a few xlOOcm" 3 in optically 
thick clouds, simila r to the average density w ithin GMCs 
in the Milky Way ijScoville fc Sanders! [l987h . The criti- 
cal density is even higher, ~ 3000 cm -3 , in optically thin 
regions. Thus, CO emission should not be detected if 
the density is lower than the densities of GMCs. Sec- 
ondly, CO molecules rapidly dissociate in the diffuse in- 
terstellar radiation field if they arc unshielded by a suf- 
ficient column of dust. It is often discussed that a visual 
exti nction of only Ay ~ 1 mag i s sufficient for shield- 
ing (|Pringle. Allen fcLu bow 200l|); however, this is true 
only at high densities (10 3 cm -3 ; Ivan Dishoeck fc Blackl 
1988). A higher Ay is necessary at lower densities, since 
the rate of molecular formation must be rapid to counter- 
balance the rapid dissociation. At the typically densities 
of a few 10 2 cm -3 expected for molecular gas based on 
Galactic studies, several mag of visual extinction are re- 
quired and the CO emitting gas must reside in GMC-like 
structures during passage trough the inter-arm areas. 

Lastly, we address the nature of the GMAs: are they 
distinct clouds, and simply an extension of GMC mass 
spectrum, or the confusion of many GMCs concentrated 
by orbit crowding in the spiral potential but unresolved 
due to the limited spatial resolution. The gas surface 
densities within the GMAs is 200-1000M Q pc~ 2 , gener- 
ally higher than that in Galactic GMCs (~ 170M Q pc" 2 ; 
iSo'lomon et ail fl987f) . The FWHM thickness of the 
molecular gas disk in M51 is unknown b ut in the Galaxy 
it is - 120 pc (jScoville fc Sanders|[l987l) . so the average 
gas densities within the G MAs (40-200 cm" 3 ) ar e similar 
to typical GMC densities ([Solomon et al.lll987t ). Thus, 
the molecular gas would continuously fill the entire vol- 
ume within GMAs. We conclude that the GMAs are 
most likely not just confusion of multiple GMCs, but 
instead must be discrete structures. The GMCs coag- 



ulate and form GMAs in spiral arms. Such mass con- 
centrations would be susceptible to kinematic fragmen- 
tation in high shear gradients across the spiral arms. In 
fact, the shear gradients are 200-600 km s _1 kpc" 1 (Fig- 
ure [TJd) , large enough to pull apart the GMAs against 
their self-gravity; the shear timescale, defined as the in- 
verse of Oorts A-constant, is comparable to the free- 
fall timescale (~ 5-10 Myr). In addition the GMAs in 
our images are in virtually all instances elongated along 
the spiral arms, instead of the round shapes expected if 
they are strongly self-gravitating. Formation of massive 
GMAs without an aid of gravity is seen in theoretical 
mode ls due to GMC agglome ration in the spiral density 
wave (jDobbs fc Bon ncll 2006) , and due to hydrodynamic 
instabilities triggered by strong shear upon e ntering the 
spiral arms (jWada fc Kod a 2004; Wada 20(1). 

The remnants of fragmented GMAs are seen in the 
inter-arm regions. Optical and infrared images show 
spur structures e merging from the spiral arms into the 
inter-arm regions (|La Vigne. Vogel fc Ostriken 120061 ) as 
dark filamentary lanes originating on the outside (down- 
stream) of spiral arms (Figure [ 3. A few spu rs have also 
been detected in CO line (|Corder et al.ll20(M) and in Fig- 
ure [3] here. They are formed by fragmentation of GMAs 
leaving the arms into the inter-arm regions. With the ob- 
served shear motions, GMAs would naturally shear into 
such filamentary spur structures. The total masses of 
spurs are a few xlO 6_7 M0, approaching the masses of 
typical GMAs (~ 10 7 M Q ). 

These new observations suggest that the evolution of 
the dense ISM in M51 is dynamically driven. The GMCs 
coalesce into GMAs as they flow into the spiral arms and 
the cloud orbits converge in the spiral arm. The very 
massive GMA seen in the spiral arms must have life- 
times comparable to the time needed to cross the arms 
(~ 2 - 5 x 10 7 yrs). On leaving the spiral arms, these 
GMAs are fragmented by the strong shear motions and 
then ejected into the downstream inter-arm regions as 
lower mass GMCs. Our new observations reveal over a 
hundred of the most massive GMCs (> 4 x 1O 5 M ) in the 
interarm regions. These GMCs are detected throughout 
the interm areas and therefore must have lifetimes com- 
parable with the interarm crossing time of 10 s yrs. The 
majority of the interarm molecular gas is not resolved 
at our current detection limit ~ 4 x 10 5 M Q - thus the 
fragmentation of GMAs must proceed to even lower mass 
GMCs. It is clear that most of the molecular gas is con- 
verted not simply to HI but to GMCs, since the molecules 
dominate the overall gas abundance in the major part of 
the disk (central ~ 12 kpc). The lifetimes of the lower 
mass GMCs require more sensitive observations. 
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Fig. 1. — (a) Integrated intensity map of CO(J=1-0) emission of the entire disk of M51. The 6.0' X 8.4' region was mosaiced in 151 
pointings at 4 resolution with the CARMA interferometer. The total power and short-spacing data are obtained with the On-The-Fly 
mapping mode of the BEARS multi-beam receiver on the Nobeyama Radio Observatory 45m telescope (NR045). The CARMA and NR045 
data are combined in the Fourier space. The map clearly detect GMCs over the entire disk for the first time, including both the prominent 
spiral arms and inter-arm regions, (b) Velocity field. Significant shear motions are seen at tangential positions (PA of the disk kinematic 
major axis is -lldeg) 



6 



Koda et al. 




Fig. 2. — (a) Distribution of GMCs (10 5 ~ 6 Mp) and giant molecular asso ciations (GMAs; > 10 7 M Q ) in M51. The GMCs/GMAs are 
identified with the CLUMPFIND algorithm (Williams, dc Gcus & Blitz 1994), down to 4<r significance, corresponding to the typical GMC 
mass in the Milky Way (4 X 10 5 Mq; Scovillc & Sanders 1987). The small green circles includes only the GMCs with mass above 4 X 10 5 Mq . 
The GMAs are seen only in the spiral arms, suggesting that they are assembled and broken up as the gas flows through the spiral arms. 
Numerous GMCs are still seen in the inter-arm regions, indicating that they survive while crossing the inter-arm regions, (b) Nobeyama 45 
m telescope CO(J=1-0) map. The circle at the lower-left corner is the 22" beam. The white contour around the emission indi cates roughly 
the coverage of CARMA observations, (c) Molecular gas fraction defined as 2n# »/( n H +2nH 2 )> calculated with Hj data from lBraun et aLl 
(2007). The fraction is high within the major part of the disk and does not change azimuthally, indicating that the gas stays molecular 
over a revolution. 




Fig. 3. — Spurs in inter-arm regions. Spurs that originate from the spiral arms and extend into the inter-arm regions are seen as optical 
extinction in the B-band image from the Hubble Space Telescope (HST) archive (color). Molecular gas (contours) traces the spurs, possibly 
the fragmented remnants of GMAs. Contours are at 3, 5, 9, 13<r in each velocity channel, and the lowest contour is presented with a thick 
line. The spurs have very narrow line widths and are seen clearly in channel maps, but not as much in an integrated intensity map. We 
therefore overlayed the contours of all channels with emission exceeding 3 a. 



